Magnetic instability in a dilute circular rarefaction wave 
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The generation of a magnetic field in a circular rarefaction wave is examined in form of a 2D 
particle-in-cell (PIC) simulation. Electrons with a temperature of 32 keV are uniformly distributed 
within a cloud with a radius of 14.2 electron skin depths. They expand under their thermal pressure 
and carry with them the cold protons, which are initially concentrated in a hollow ring at the 
boundary of the electron cloud. The interior of the ring contains an immobile positive charge 
background that compensates for the electron charge. The protons expand in form of a circularly 
symmetric rarefaction wave and they extract energy from the electrons. A thermal anisotropy of the 
electrons develops and triggers through a Weibel-type instability the growth of TM waves within the 
plasma cloud, which acts as a waveguide. The changing cross section of this waveguide introduces 
a coupling between the TM wave and a TE wave and in-plane magnetic fields grow. The relevance 
of the simulation results to a previous experimental study of a laser-ablated wire is discussed. 

PACS numbers: 52.65.Rr, 52.35. Qz, 52.38.Fz 



I. INTRODUCTION 

Magnetic instabilities in expanding plasmas or in den- 
sity gradients are of interest for laser fusion, where mag- 
netic fields can reduce the particle's mobility [1], and 
for astrophysical plasmas, where they can generate mag- 
netic fields from noise levels. Magnetic field generation in 
the interstellar medium (ISM) Q by the interplay of the 
electrons of the hot ionized medium with the spatially 
nonuniform (ISM) plasma or the amplification of mag- 
netic fields in the turbulent plasma close to supernova 
remnant (SNR) shocks are examples. 

The growth of magnetic fields in rarefaction waves on 
an electron timescale has recently been observed 

experimentally [Toj and with particle-in-cell (PIC) sim- 
ulations 0, El- The instability is driven by a thermal 
anisotropy and is thus similar to the Weibel instabil- 
ity [Ea - tl6| |. It is the result of the electron's slowdown 
by the ambipolar electrostatic field, which is sustained 
by the plasma density gradient of the rarefaction wave. 
This electrostatic field counteracts the charge separation 
along the plasma density gradient that arises from the 
difference in the thermal speeds of electrons and ions. It 
accelerates the ions and slows down the electrons along 
this direction, which generates the thermal anisotropy. 

Previous simulations [(1 0,112 have considered systems, 
in which hot electrons accelerate an equal number of ini- 
tially cold ions. Here the acceleration of protons, which 
are distributed in form of a hollow ring, by the ambipo- 
lar electrostatic field is examined with a PIC simulation. 
This proton ring distribution is a good approximation of 
the cross section of the ions of a laser-heated wire, as long 
as it is located far away from the laser impact point. 
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The ablation of the wire and the resulting magnetic 
instabilities have been examined experimentally in Ref. 
[To| . The wire in the experiment is composed of heavy 
ions, which can not be accelerated to high speeds by the 
expanding electrons. The light protons that feed the rar- 
efaction wave originate from surface impurities and they 
are approximated here by the hollow ring distribution. 
The reduced number of mobile protons dilutes the rar- 
efaction wave. This is also observed in theexperiment 
where its number density is ~ 10 18 cm~ 3 [10(, which is 
well below the solid ion number density of the wire. How- 
ever, only the proton distribution in our simulation is re- 
alistic with respect to the experiment. The electrons have 
the temperature 32 keV and are spatially uniform within 
the plasma cloud, while the electrons in the experiment 
have MeV energies and are confined to the wire's surface 
[T?| . Reducing the electron temperature and distribut- 
ing the electrons over a wide interval is computationally 
efficient and it ensures that the overall thermal energies 
in the experiment and simulation are comparable. 

The purpose of our simulation is threefold. We want 
to determine if the gradient- driven magnetic instability 
always develops in the rarefaction wave or if competing 
instabilities can outrun it. Secondly, we want to deter- 
mine if the lower proton density results in a weaker ther- 
mal anisotropy of the electrons and, thirdly, if and how 
the expansion of the dilute rarefaction wave differs from 
the dense one in Ref. pd| . 

Our results are as follows. The protons at the front of 
the dilute rarefaction wave are accelerated to about the 
same speed within the same time interval as those in the 
dense one in Ref. [ll]. This confirms our expectation. 
The electrostatic potential of the plasma with respect to 
the surrounding vacuum is fully determined by the ther- 
mal pressure of the electrons, which is the same here and 
in our simulation in Ref. The stronger reduction 

of the electron's thermal energy in Ref. compared 
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to the one here does not lead to detectable differences in 
the plasma expansion. The same magnitude of the elec- 
trostatic potential and the therefrom resulting equal de- 
celeration of the electrons in the direction of the density 
gradient imply that the electron's thermal anisotropics 
here and in Ref. [ll[ are equally strong. However, the 
much lower electron number density in the rarefaction 
wave we model here delays the onset of the gradient- 
driven magnetic instability. It develops instead in the 
dense core of the plasma cloud and the strong magnetic 
fields diffuse out into the rarefaction wave. 

The magnetic fields grow in the radial interval of the 
plasma cloud in which the temperature anisotropy and 
the number density of the electrons are large. The mag- 
netic field source is thus a Weibel-type instability and 
not the competing thermoelectric instability 18]. The 
latter is inefficient in our case study, because the electron 
density gradient is aligned with the electron temperature 
gradient and because the expanding rarefaction wave is 
circularly symmetric. 

We also observe here the same secondary magnetic 
instability as in Ref. [ll|, which yields the growth of 
in-plane magnetic fields. Our present simulation setup 
confines the secondary wave in the plasma cloud's core, 
while it developed in the expanding rarefaction wave in 
Ref. [TTj . This confinement simplifies the interpretation 
of the data. Our results suggest that the in-plane mag- 
netic fields are generated by a mode conversion of a TM 
wave into a TE wave. This is a well-known process in 
waveguides with a variable cross section, which are im- 
portant in antenna theory [2l| . The variable cross section 
is here the result of the proton expansion. 

The structure of this paper is as follows. The PIC sim- 
ulation scheme and the initial conditions are summarized 
in Section 2. Section 3 presents the simulation results, 
which are discussed in Section 4. 



II. THE SIMULATION CODE AND THE 
INITIAL CONDITIONS 

A particle-in-cell (PIC) code approximates a plasma 
by an ensemble of computational particles (CPs) and it 
uses their collective charge distribution p(x) and current 
distribution J(x) to evolve in time the electromagnetic 
fields on a spatial grid. The electric E and magnetic B 
fields update in turn the momentum of each CP through 
the relativistic Lorentz force equation. The PIC scheme 
is discussed in detail in Ref. (19| . 

Most codes evolve the electromagnetic fields through 
the discretized forms of the Ampere's and Faraday's laws 
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is used to update the momentum pj of the j parti- 
cle of species i. The collective behavior of the ensem- 
ble of the CPs of species i approximates well that of a 
plasma species, provided that the charge-to-mass ratio of 
the plasma particles equals the ratio qi/m.i of the CPs, 
that the plasma is collision-less and that the statistical 
representation of the computational plasma is adequate. 
We use the numerical scheme discussed in Ref. [20]. 

Our initial conditions and their connection to the ex- 
periment are visualized in Fig. Q] The axis of the wire on 
the left hand side is parallel to z. The simulation resolves 
a cross-section of this wire in the x-y plane with an origin 
x = 0, y = in the center of the wire. This cross-section 
has a z-coordinatc that is sufficiently far away from that 
of the laser impact point so that we do not have to model 
the laser pulse in the simulation. 

In the experiment, the hot electrons stream uniformly 
from the laser impact point along the wire's surface [17| 
to the cross section that corresponds to our simulation 
plane. We thus approximate the wire's cross section by 
the circular plasma cloud shown on the right hand side 
of Fig. Q] Hot electrons fill the entire cross section of 
the plasma cloud with radius r = rw- Their spatially 
uniform number density n$ within the cloud gives the 
plasma frequency u> p . The mobile protons fill a hollow 
ring with the outer radius r = rw and with the inner ra- 
dius r — 0.95rw The interior r < 0.95rw of the hollow 
ring contains an immobile positive charge background. 
The charge density of the electrons equals at any loca- 
tion x, y with x 2 + y 2 < r 2 that of the positive charge 
carriers and the mean speeds of both mobile species are 
zero. No net charge and no net current are thus initially 
present. All initial electromagnetic fields are thus set to 
zero. The cloud is immersed in a vacuum and the bound- 
ary conditions of the simulation box are periodic. 

The outer cloud radius rw = 14.2A e , where A e = c/ui p 
is the electron skin depth within the cloud. The temper- 
ature of the relativistic Maxwellian distribution that rep- 
resents the electrons is 32 keV and their thermal speed is 
v e = c/4. The bulk of the electrons thus moves at nonrel- 
ativistic speeds, which allows us to easily decompose their 
radial and azimuthal velocity components. The electron 
Debye length A^ = v e /uj p equals 5.6A,,;, where A x is the 
side length of a quadratic grid cell. We use the correct 
mass ratio between the electrons and the protons in the 
ring distribution. The proton temperature is 10 eV and 
their thermal speed v p = 3.1 x 10 4 m/s. 

The quadratic area L x L of the 2D simulation box with 
the side length L = 75.5A e is resolved by N g — 1700 grid 
cells in the x and y directions. We represent the electrons 
by N e = 8 x 10 8 CPs and the protons by N p = N e 
particles. Since the protons occupy a smaller volume, 
their numerical weight is lower. A time interval tu> p = 800 
is resolved. We normalize time to lj p , space to the initial 
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FIG. 1: The initial conditions: The experimental setup is 
shown to the left, where we assume that a long wire is aligned 
with z. The initial electron distribution in the simulation is 
shown to the upper right and the proton ring distribution 
(dark gray shade) and the positively charged immobile back- 
ground (light gray shade) to the lower right. 



cloud radius r\y and the electron and proton velocities 
to their respective initial thermal speeds v e and v p . The 
electric and magnetic fields are normalized to m e oj p c/e 
and m e ujp/e, respectively. 



III. SIMULATION RESULTS 

The kinetic energy of electrons with the mass mi is 
K\{t) — rriic 2 Sfci P\f — 1]) where the summation is 
over all computational electrons with the Lorentz fac- 
tors Tj and Kq = Ki(0). The kinetic energy of the 
mobile protons with mass iri2 = 1836 mi is K2(t) — 
m 2C 2 J2j=i [Fj — 1]. The energy of the in-plane elec- 
tric field is E E± (t) = (Al/2e )J2?] =1 E 2 (i.j.t), where 

E p (i,j,t) = [E%(i,j,t)+E*(i,j,t)] 1/3 . The energy 
density of the out-of-plane magnetic field Esz{t) = 

(A* 2 Mo )E3 : niu.j.0. 

Figure [5] shows their time evolution. The electrons sus- 
tain a rapid energy loss during < t < 20. This energy 
is transfered to the ambipolar electric field, which grows 
and saturates during this time. This field accelerates the 
protons and the electrons have transfered about 10% of 
their initial energy to them at t = 800. The initial os- 
cillations of Ee± have damped out at t ~ 200 and Ee± 
reaches a steady state value of ~ 10 -2 -K"o- Eb z grows ini- 
tially slowly. The faster growth of the magnetic energy in 
the time interval 400 < t < 700 is followed by its satura- 
tion. The magnetic energy remains well below that found 
in Ref. [ll| and about two orders of magnitude below the 
electric one. We will examine now in more detail the field 
and particle distributions at the time t = 27 when the 
electrostatic field reaches its peak value, at t = 500 when 
the magnetic field grows fastest and at t = 800 when the 
magnetic field saturates. 

The modulus of the in-plane electric field at t = 27 
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FIG. 2: Energy densities in units of the initial electron ther- 
mal energy: The electron energy (a), the proton energy (b), 
the energy density of the in-plane electric field (c) and of the 
out-of-plane magnetic field (d). 



is shown in Fig. [3J The electric field has the expected 
circular symmetry. It peaks at r » rw It gradually de- 
creases for increasing values of r and reaches noise levels 
at r « 2rw- Circular electric field oscillations are visible 
in the cloud's core r < rw- Since the positive charged 
background is immobile in this region, these electrostatic 
waves must be Langmuir waves. The magnetic field re- 
mains at noise levels at this time (not shown). 

The phase space density distributions f e (r,v r ) and 
f p {r,v r ) of electrons and protons, which are functions 

of the radius r — (x 2 +y 2 ) 1 ^ 2 and of the radial veloc- 
ity component v r — (v 2 +v 2 ) 1 ^ 2 , are displayed in Fig. 
21 They are derived as follows. The circular symmetry 
of the cloud and the energy exchange between electrons 
(s=e) and the mobile protons (s=p), which is at least 
initially limited to the x,y-plane and a function of the ra- 
dius, imply that bi-Maxwellian distributions f s (r, v r ,v z ) 
will develop for both species. We neglect the v z direction 



FIG. 3: (Color online) The 10-logarithmic modulus of the 
electric field \E p (x,y)\ sampled at the time t = 27 (Enhanced 
online). 

and define the average over the azimuthal angle p 

f s (r,v r ) = (2irr)- 1 [ f s {r,v r )rd Pl (4) 

which greatly improves the visualized dynamical range 
of the phase space densities. We normalize f s (r,v r ) by 
its maximum value at t = 0. The resulting phase space 
densities are thus constant as a function of r at t = 
within the radial interval the mobile species occupy. 

The protons have been accelerated to several times v p 
at t = 27 and the proton mean speed increases with the 
radius. The protons have, however, not moved far beyond 
rw during this short time. The fastest electrons have 
reached a radius r w 2.3rw- Their density profile shows 
a straight line from r w 1.5rw and v r « to r ~ 2.3rw 
and v r ~ 4. These are the electrons that escaped from 
the cloud before the ambipolar electric field has fully de- 
veloped. The linear density profile simply reflects that 
faster electrons have propagated to larger radii during 
t = 27. The ambipolar electric field, which has been 
built up after the first electrons escaped into the vacuum 
(See Fig. (He)), affects the electrons at lower r, and we 
can observe a decrease of the peak electron speed as we go 
from r 2.3rw to r rs l.brw- The electrons lose kinetic 
energy as they overcome the electrostatic ambipolar po- 
tential. This potential is also responsible for the drastic 
drop of the electron number density at r w rw ■ It is evi- 
dent from Fig. 0]that the electron charge for r > 1.03r\y 
is not compensated by a proton charge. The electric field 
at r > 1.03rw in Fig. f3]is thus sustained by the electron 
sheath. 

The density distributions of the electrons and the mo- 
bile protons are shown in Fig. [SJ They are obtained 
from the integration of the phase space densities in 
Fig. 0] over v r . They are normalized to their initial 



FIG. 4: (Color online) The 10-logarithmic phase space densi- 
ties f r (r,v r ) of the electrons (a) and f r (r,v r ) of the protons 
(b). The densities are normalized to their initial value. The 
simulation time is t — 27 (Enhanced online). 




FIG. 5: The electron density n e (r) and the proton density 
rip(r) (dashed curve) sampled at the time t = 27 and normal- 
ized to the respective initial density values. 



value. The density distributions confirm that an elec- 
tron sheath has formed outside the expanding proton 
cloud. The proton density has spread from the initial 
interval 0.95 < r/r w < 1 to 0.95 < r/r w < 1.03. Note 
that we have not plotted here the charge density inside 
the proton ring distribution, The density equals 1 in our 
normalization for r < 0.95rw and it ensures an overall 
charge neutrality of the plasma. The density distribu- 
tions of the mobile species demonstrate that the interval 
0.98 < r/rw < 1-02 is positively charged. 

The electron and proton phase space density distribu- 
tions / e (r, v r ) and f p (r, v r ) at the time t = 500 are shown 
in Fig. [BJ The majority of the electrons in Fig. [Ha) is 
confined by the immobile positive charge background in 
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FIG. 6: (Color online) The lOTogarithmic phase space densi- 
ties fr(r,Vr) of the electrons (a) and fr(r,v r ) of the protons 
(b). The densities are normalized to their initial value. The 
simulation time is t = 500. 
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FIG. 7: The electron density n e (r) and the proton density 
n p (r) (dashed curve) sampled at the time t = 500 and nor- 
malized to the respective initial density values. 




the interval r < 0.95rw- Their phase space density and 
the characteristic electron speed decrease drastically at 
r ps rw but remain relatively high up to r ~ l.brw- The 
phase space density decreases by another two orders of 
magnitude as we go to even larger radii. The radial in- 
terval 1 < r/rw < 1-5 with the elevated electron phase 
space density and mean speed coincides with the radial 
interval that is occupied by the distribution of mobile 
protons. An almost closed circular proton phase space 
structure is present at r ps 1.03rw in Fig. EJb), which is 
similar to an ion phase space hole, and it must be related 
to a local excess of negative charge. The profile of the 
proton phase space density distribution in Fig. &b) in- 
creases linearly with the radius for r > 1.05r\y, which is 
characteristic for a rarefaction wave. The protons reach 
a peak speed « 200v p , which is comparable to v e /10 and 



about the same as that in Ref. [11] . Only a small frac- 
tion of the protons reaches this speed, which limits the 
loss of electron thermal energy at this time (See Fig. [5]). 
The low electron phase space density for r > 2rw implies 
that all electron processes close to the boundary are slow 
and do not carry much energy, which justifies our choice 
of periodic boundary conditions. 

The densities of the mobile particle species at t = 500 
are shown in Fig. [71 The electron density decreases by an 
order of magnitude close to the boundary of the immobile 
positive charge background at r = 0.95rvp. It increases 
again for r > rw and reaches a local maximum at around 
r ps 1.15rw, close to the peak of the mobile proton's 
density. Both densities decrease gradually beyond this 
radius. The proton density falls off steeply at its front at 
r l.5r w (See also Fig. E^b)). 

Figure [5] shows the in-plane electric field distribution 
at t = 500. The electric field modulus shows a more 
complex pattern than at the earlier time. The electric 
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FIG. 8: (Color online) The 10- logarithmic modulus of the 
electric field \E p (x,y)\ sampled at the time t = 500. 



field peaks at r ~ rw and it is sustained by the electron 
density gradient close to the boundary of the immobile 
positively charged background at r = 0.95rw- The elec- 
tric field modulus goes through a minimum at a slightly 
larger radius. The reason for this radial oscillation is that 
we have two regions with an excess of positive charge, 
which are separated by a radial interval with a negative 
excess charge at r « rw (See Fig. [JJ. An electric field 
with a significant amplitude modulus that extends over a 
large radial interval is observed at 1.1 < r/rw < 2. The 
electric field modulus in this band peaks at r ps l.5rw, 
which coincides with the front of the proton distribution 
in Fig. [7J It thus corresponds to the ambipolar electric 
field driven by this density gradient, and it decreases for 
larger radii where we only find electrons. The electric 
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field modulus reaches its minimum at r ~ 2ryy- 

Figure |H] displays the distribution of B z at t — 500. 
We observe a strong localized magnetic field structure. 
Its amplitude peaks at r « 0.8rw and it decreases as we 
go to r ss rw and beyond. The amplitude oscillates once 
in the azimuthal direction and the associated wavelength 
A p ~ 2irrw is thus much larger than the one in the radial 
direction, which we can estimate as follows. We find a 
magnetic field maximum at x ~ — O.Srw and y ~ O.Srw 
and a minimum at x ~ —0.6r\y and y ~ 0.6r^. The 
wavelength of the radial oscillation is thus A r sa 0.5r\y- 
The magnetic amplitude outside this radial interval is at 
noise levels. The magnetic noise is distributed over the 
entire simulation box, while the magnetic structure is lo- 
calized in a small radial interval. This explains why we do 
not observe a more pronounced growth of the magnetic 
field energy in Fig. [UJd). 

It is instructive to compare the radial interval, in which 
the magnetic field grows, with the one that shows a 
thermal anisotropy. We determine for this purpose the 
thermal energy densities of the electrons in the radial 
and azimuthal directions. We consider only the in-plane 
component of the speed v pj - = (v x ,v y )j of the j th com- 
putational electron. The radial component of the ther- 
mal energy K r j = t£ • is computed from the projection 
Vr,j = v p.j ' v jl T ji where Tj is the position vector of the 
electron in circular coordinates, and K p j = vz j — K r j. 
The partial thermal energies and the anisotropy A are 
then obtained from the summations 

K r (i5 r ) = K r jSjj (5) 

K p {i8 r ) = Y i K p , j 6 id (6) 

A = K r /K p , (7) 

where Sij = 1, if rj falls into the i th radial bin of width 
S r , and zero otherwise. We thus obtain a histogram A(i) 
of the radial distribution of the thermal anisotropy. The 
anisotropy A(r) is compared with the magnetic energy 

P Bz (nS r ) = (A3/2 Mo ) ]T B\{i - N g /2,j - N g /2)S Lj>n , 

(8) 

where Sij, n — 1 if I(i 2 + j 2 ) — n 2 , with I being a round- 
off operation and = 5 r . This azimuthal integration, 
rather than the azimuthal average, emphasizes magnetic 
fields at larger radii. 

Figure [TOl demonstrates that the magnetic field starts 
to grow at t w 200, when an anisotropy A(r) < 1 has 
formed that is sufficiently strong and wide. The magnetic 
field grows initially in the interval 0.7 < r/rw < 1-1, but 
it expands later on in both radial directions. Its front 
reaches r ~ 1.5rw at t = 500, which coincides with the 
tip of the proton distribution in Fig. [7] It is thus confined 
to within the rarefaction wave. The radial interval where 
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FIG. 10: (Color online) The ratio A between the mean ra- 
dial energy and the mean perpendicular energy of electrons 
is shown on a linear color scale in (a). The spatio-temporal 
evolution of the magnetic energy density is shown on a 10- 
logarithmic color scale in (b). 



Pbz ~ 1 is initially stationary, but the magnetic field 
distribution changes after t « 500. What appears to be 
a sidelobe develops at r k lArw- This change takes 
place on a few tens of lo~ 1 . The plasma frequency in the 
rarefaction wave is about Wp/4 at t = 500 (See Fig. [7]) 
and the growth time of the sidelobe is thus faster than 
that expected from any instability. Indeed, the on-line 
enhancement of Fig. |H] shows that the magnetic i? 2 -field 
leaks out from within the cloud into the rarefaction wave. 

Figure [Til shows its distribution at t = 800. The ex- 
trema of B z are located within r — rw and the amplitude 
of B z is constant as a function of r for 1 < r/rw < 1-5. 
The magnetic energy has been confined to within r sa ry/ 




FIG. 11: (Color online) The magnetic field amplitude B z (x, y) 
sampled at the time t — 800. Overplotted is a circle of radius 
r = r W - 




FIG. 12: (Color online) The magnetic field amplitude 
\B p (x,y)\ sampled at the time t = 800. Overplotted is a 
circle of radius r — rw ■ 



until t sa 600 and it forms a TM mode until then (See 
Fig. [9]). The rapid expansion of the magnetic energy af- 
ter t « 600 implies that this TM mode suddenly expands. 

We can understand the circular plasma cloud as the 
cross section of a cylindrical waveguide with an axis that 
is aligned with z. The expansion of the mobile protons 
imply that the radius of the cross-section of this waveg- 
uide increases in time. It is well-known that changes in 
the radius of a waveguide induce a coupling of TM and 
TE modes [2lJ. The magnetic field of a TE mode would 
be oriented in the simulation plane. 

Figure [12] evidences that a TE wave is indeed present 
at t = 800. A structure is visible in the distribution of the 

1 /2 

in-plane magnetic field B p = (B^, + By) . It is located 
in the same radial interval as the magnetic field of the TM 
mode. The magnetic field patterns in B z and B p have the 
same azimuthal wave number and both are phase-shifted 
by 90° relative to each other. The magnetic energy of 
the TE mode grows after t = 800 to values exceeding 
that of the TM wave (not shown) . We do not discuss the 
evolution for t > 800. The TE wave expands out to the 
boundaries shortly after this time, which results in finite 
box effects, while the TM wave remains confined by the 
rarefaction wave as discussed previously [ll[ . 

The different behaviour of the TM and TE modes can 
be explained by the different plasma response to their 
electric fields in our 2D box geometry. An electric field 
orthogonal to the simulation box plane can not result in 
charge density modulations, since we do not resolve the 
z-direction. It thus only affects the current distribution. 
An in-plane electric field does, however, modulate also 
the charge density. The dilute electron plasma between 
the front of the rarefaction wave and the boundaries can 
not support strong charge density waves, but it can easily 
support the large currents from the TE wave and the 
latter can expand more easily. 



IV. DISCUSSION 

We have modeled here with a 2D PIC simulation the 
expansion of a circular plasma cloud into a vacuum, 
which has been driven by the thermal pressure of the 
electrons. It is a follow-up study of a previous simulation 
experiment. It aimed at explaining the growth of mag- 
netic fields in the rarefaction wave, which is g enerated 
by the ablation of a wire by a laser pulse [ltj | ■ It con- 
sidered the expansion of a circular plasma cloud, which 
consisted of spatially uniform hot (32 keV) electrons and 
cool (10 eV) protons. Here we have confined the mobile 
protons to the border of the plasma cloud. This hol- 
low ring distribution is a more accurate approximation 
of the experimental conditions. The rarefaction wave ob- 
served in Rcf. flo| contains primarily the light ions from 
the surface impurities, which have been ionized by the 
strong surface electric field and current [17]. However, 
computational constraints require us to represent here 
the electrons as a hot (32 keV) species that is uniformly 
distributed over the entire plasma cloud. The electrons 
in the experiment reach MeV temperatures, but they are 
confined to the wire's surface. Choosing cooler electrons 
reduces the difference between the electron and proton 
Debye lengths, which is computationally efficient, while 
it ensures that the thermal energy that drives the expan- 
sion is comparable in simulation and experiment. 

Our results are as follows. The ambipolar electric field 
driven by the electron's thermal expansion results in the 
formation of a rarefaction wave. The fastest protons 
reach a speed that is comparable to about a tenth of the 
electron thermal speed, which equals the value observed 
in Ref . [ll[ . The proton acceleration is thus not affected 
by the choice of the initial proton distribution. However, 
the density of the rarefaction wave is limited by the num- 
ber of available mobile protons. The spatially uniform 
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proton distribution in Ref. [ll| provided a continuous 
feed of mobile protons, while the number of mobile pro- 
tons we introduce here is limited. The plasma density in 
the rarefaction wave we observe here is thus lower than 
that in Ref. [ll|. The therefrom resulting lower plasma 
frequency implies a slowdown of the instabilities in the 
rarefaction wave. 

Although the thermal anisotropy in the electron dis- 
tribution here and in Ref. [ll| has been comparable, the 
Weibel-type instability could not develop here during the 
simulation time. The instability developed instead in the 
dense core of the plasma cloud and the magnetic field 
diffused out into the rarefaction wave. This magnetic 
diffusion is not compatible with the magnetic field struc- 
tures observed in Ref. [13], which can be explained bet- 
ter by an instability in the rarefaction wave. Remarkably 
the more realistic proton distribution we use here yields 
results that are less representative for the experimental 
ones than those of the simulation in Ref. [I(|. This is 
probably a consequence of the spatially uniform electron 
distribution. The MeV electrons form a thin layer on 
the wire's surface in the experiment and the magnetic 
instability can not form in it. It can only form in the 
rarefaction wave, which maintains a significant thermal 
anisotropy over a large radial interval. Future studies 
thus have to confine the electrons to a smaller radial in- 
terval. 



The present simulation sheds light on the mechanism 
by which the in-plane magnetic fields grew in Ref. [Tl| . 
The plasma in the radial interval r < 0.95r\y with the 
immobile positive charge background is equivalent to a 
cylindrical waveguide and the thermal anisotropy of the 
electrons results in the growth of a TM wave inside the 
wave guide. The plasma of the dilute rarefaction wave 
is a perturbation of the waveguide's cross section. A 
varying radius of a waveguide can couple TM and TE 
modes [2lT | . Here this coupling results in the growth of 
in-plane magnetic fields. A TE mode can probably not 
be driven by a plasma instability. The electromagnetic 
forces and the plasma flow are confined to within the 
simulation plane. No current can thus develop in the 
orthogonal direction. Orthogonal plasma currents are, 
however, needed to maintain the in-plane magnetic field. 
A wave coupling between a TM and a TE wave drives 
orthogonal electric fields, which can accelerate electrons 
in this direction. Future work will address this wave cou- 
pling with a larger simulation box. The evolution of the 
TE wave can then be examined for a longer time without 
finite box effects. 
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